Cross sections for hard exclusive electroproduction of ir + mesons on a hydrogen target 
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The exclusive electroproduction of 7r + mesons was studied with the Hermes spectrometer at the 
Desy laboratory by scattering 27.6 GeV positron and electron beams off an internal hydrogen gas 
target. The virtual-photon cross sections were measured as a function of the Mandelstam variable t 
and the squared four momentum — Q 2 of the exchanged virtual photon. A model calculation based 
on Generalized Parton Distributions is in fair agreement with the data at low values of \t\ if power 
corrections are included. A model calculation based on the Regge formalism gives a good description 
of the magnitude and the t and Q 2 dependences of the cross section. 
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PACS numbers: 13.60.-r, 13.60.Lc, 13.85.Lg, 14.20.Dh, 14.40.Aq 



The interest in hard exclusive processes has grown 
since a Quantum Chromodynamics (QCD) factorization 
theorem was proven for the hard electroproduction of 
mesons by longitudinal photons 0]. It was shown that 
at leading order in the strong-coupling constant as and 
in the fine-structure constant a the amplitude for such 
reactions can be factorized into three parts (see left panel 
of Fig. [T]): a hard lepton-scattering part, which can be 
calculated in Quantum Electrodynamics and in perturba- 
tive QCD (pQCD), and two soft parts that parametrize 
the structure of the target nucleon by Generalized Par- 
ton Distributions (GPDs) @, 0, Q and the structure of 
the produced meson by a distribution amplitude. The 
GPDs offer the possibility to reveal a three-dimensional 
representation of the structure of hadrons at the partonic 
level, correlating the longitudinal momentum fraction to 
transverse spatial coordinates [ 1 , B- For recent 

theoretical reviews, see Refs. [I(llllill2j. 

The amplitude for exclusive electroproduction of 
mesons with specific quantum numbers is described by 
a particular combination of GPDs. At leading twist, ex- 
clusive vector-meson production is sensitive to only un- 
polarizcd GPDs (H and E), while pseudoscalar-meson 
production is sensitive to polarized GPDs (H and E) 
without the need for a polarized target or beam. More- 
over, for 7r + production, the pseudoscalar contribution 
involving E dominates at small momentum transfer to 
the target as it contains the i— channel pion-pole contri- 
bution. The complete amplitude of that contribution to 
7r + production also contains the pion charge form fac- 
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At moderate virtuality of the exchanged photon, ncxt- 
to-leading-order (NLO) in as and higher-twist correc- 
tions to the pQCD leading-order amplitude can con- 
tribute. Two different types of higher-twist corrections, 
jointly denoted by the term "power corrections" , have 



been estimated in the case of 7r + production [151 ]: one 



arising from the intrinsic transverse momentum of the 
partons, and the other resulting from the soft-overlap 
diagram. In contrast to the leading-order perturbative 
mechanism, the latter, shown on the right panel of Fig.JTJ 
does not proceed through one-gluon exchange. Although 
significant NLO corrections have been calculated [l6l.[l7l]. 
the higher-twist corrections dominate for the virtuality of 
the exchanged photon of the present data. 

Alternatively, exclusive processes can be described at 
the hadronic level within the Regge formalism (e.g. [l8j|) 
where the interaction between the virtual photon and 
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FIG. 1: Leading-order (left panel) and soft-overlap (right 
panel) diagram for the exclusive 7r + electroproduction ampli- 
tude. 



the nucleon is described in terms of Regge-trajectory ex- 
changes in the t-channel. This approach can describe 
photoproduction [3] and electroproduction [2(j of pseu- 
doscalar mesons above the resonance region. In the latter 
case, meson production by both transverse and longitu- 
dinal photons contribute and can be calculated within 
the Regge formalism. 

This letter reports a measurement of the virtual- 
photon cross section for the hard exclusive reaction ep — > 
emr + on a hydrogen target. The relevant kinematic 
variables of the process are the squared four-momentum 
— Q 2 of the exchanged virtual photon, either the Bjorkcn 
variable ig = Q 2 /2M p v (where M p is the proton mass 
and v is the energy of the virtual photon in the target 
rest frame) or the squared invariant mass of the photon- 
nucleon system W 2 , the Mandelstam variable t, and the 
azimuthal angle (f> of the pion around the virtual pho- 
ton momentum relative to the lepton scattering plane. 
Instead of t, the quantity t' = t — to was used in the 
analysis, where —to represents the minimum value of — f 
for a given value of Q 2 and xb- 

The virtual-photon cross section has been previously 
measured above the resonance region (for W 2 > 4 GeV 2 ) 
at Cea lilll. Cornell 22 1, Desy [j^ . and more recently 
at JLab j24 125I] . The present measurement extends the 
kinematic region to higher values of W and higher values 
of -t'. 

The data were collected with the Hermes spectrom- 
eter [H during the period 1996-2005. The 27.6 GeV 
Hera electron or positron beam at Desy scattered off 
polarized and unpolarizcd proton targets. Events were 
selected in which only one lepton and one positively 
charged hadron were detected and in which no additional 
cluster was recorded by the electromagnetic calorimeter. 
The Hermes geometrical acceptance of ±170 mrad hor- 
izontally and ±(40 — 140) mrad vertically results in de- 
tected scattering angles ranging from 40 to 220 mrad. 
Leptons were distinguished from hadrons with an average 
efficiency of 98% and a hadron contamination less than 
1% by using an electromagnetic calorimeter, a transition- 
radiation detector, a preshower scintillation counter, and 
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a threshold gas Ccrcnkov counter. In 1998 the thresh- 
old gas Ccrcnkov counter was replaced by a Ring Imag- 
ing Cerenkov detector {22I]. The threshold gas Cerenkov 
counter (Ring Imaging Ccrcnkov detector) provided pion 
identification in the momentum range 4.9 GeV < p < 15 
GeV (1 GeV < p < 15 GcV). For the exclusive data 
sample, the pion momentum was required to be 7 GeV 
< p < 15 GeV. For this momentum range, the pion iden- 
tification efficiency is on average 97% (99%) and the con- 
tamination from other hadrons less than 3% (2%) for the 
threshold gas Cerenkov counter (Ring Imaging Ccrcnkov 
detector). 

The kinematic requirement Q 2 > 1 GeV 2 was imposed 
on the scattered lepton in order to select the hard scat- 
tering regime. The value of W 2 was required to be higher 
than 10 GeV 2 to avoid the low-acceptance region for 
the hadron defined by the spectrometer upper angular 
limit of 220 mrad. The resulting kinematic range is 1 
GeV 2 < Q 2 < 11 GeV 2 and 0.02 < x B < 0.55. The 
mean W 2 value of the data is 16 GeV 2 . 

As the recoiling neutron was not detected, exclu- 
sive meson production was selected by requiring that 
the squared missing mass M% of the reaction ep — > 
eir + X corresponds to the squared neutron mass. Due 
to the limited experimental resolution, the exclusive ir + 
channel cannot be separated from the neighboring chan- 
nels (defined as background channels) with final states 
such as ir + + (Ntt) and ir + + (Ntttt), as their M\ values 
can be smeared into the region corresponding to exclusive 
7r + channel. In order to subtract the background chan- 
nels, which cannot be fully described by existing Monte 
Carlo simulations, a two-step procedure was developed. 
In the first step, the yield difference between ir + and tt~ 
was used. In this yield difference, exclusive ir + events 
remain as the exclusive production of 7r~ on a hydro- 
gen target with a recoiling nucleoli in the final state is 
forbidden by charge conservation. Background events 
arising from the production of neutral vector mesons 
cancel as they contribute equally to tt + and 7r~ pro- 
duction. In the second step, the Pythia Monte Carlo 
generator [28| was used in conjunction with a special 
set of Jetset [i^] fragmentation parameters tuned to 
provide an accurate description of deep inelastic hadron 
production in the Hermes kinematic domain. The dif- 
ference 1X1^+ — Nj,.-, where N represents the yield nor- 
malized by the integrated luminosity £, is well described 
by Pythia if pion momenta larger than 7 GcV are re- 
quired. The latter constraint removes mainly background 
events for which Q 2 < 3 GcV 2 . The upper panel of Fig. [2] 
shows the M\ dependence of the normalized-yield differ- 
ence — N w - for the data and for the Pythia Monte 
Carlo simulation. The Monte Carlo sample describes 
the data for values of the squared missing mass higher 
than 2 GeV 2 , i.e., outside the region corresponding to 
exclusive 7r + production. Finally, the exclusive ir + yield 
was obtained by subtracting the normalized-yicld diffcr- 
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FIG. 2: Upper panel: squared missing-mass dependence of 
the normalized-yield difference N„.+ — N„.- for data (filled 
points) and Pythia Monte Carlo (histogram) . The error bars 
represent the statistical uncertainty. Lower panel: squared 
missing-mass dependence of the normalized-yield after back- 
ground subtraction procedure. The data (filled points) are 
compared to a Monte Carlo sample for exclusive ir + produc- 
tion (histogram) normalized to the data. The inner error bars 
represent the statistical uncertainties and the outer error bars 
represent the quadratic sum of statistical and systematic un- 
certainties. The latter originate from the background subtrac- 
tion procedure. The dashed vertical line indicates the squared 
neutron mass. 



ence of the Pythia Monte Carlo from that of the data: 
N£+ = (N^+ - N 7r -) data - (N 7r+ - N^-) PYTHIA . The lower 
panel of Fig. [2] shows the peak corresponding to exclu- 
sive 7r + production resulting from this double difference. 
The peak is centered at the squared neutron mass and 
its width of 0.67 GeV 2 is consistent with that of a Monte 
Carlo sample for exclusive tt + production normalized to 
the data. The Monte Carlo, denoted as exclusive Monte 



Carlo, is based on a GPD model [15[ and will be described 
below. As Pythia does not simulate nucleon resonance 
production, the agreement of the exclusive Monte Carlo 
sample with the data for the double difference is an in- 
dication that there is very little contribution from reso- 
nant channels (n + + A for ir + and n~ + A ++ for 7r~) 
to the normalized-yield difference 1X1^+ — N„.-. In order 
to estimate the systematic uncertainty of the background 
subtraction, the Pythia Monte Carlo normalized-yield 
difference (N 7r + — N 7r -) PYTHIA was changed by the dis- 
crepancy between that normalized-yield and the data in 
the region 3 < M\ < 7 GcV 2 . The discrepancy amounts 
to between 20% and 50%, depending on the specific kine- 
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matic bin in Q 2 , xb, or t'. The largest discrepancies cor- 
respond to the least populated bins (Q 2 > 4 GcV 2 and 
-t' > 0.3 GeV 2 ). An upper limit on M 2 X of 1.2 GcV 2 
was chosen in order to optimize the combined statistical 
and systematic uncertainties for the number of exclusive 
events. The resulting relative systematic uncertainty for 
the number of exclusive events ranges from 5% to 20%. 
Within the M§- limit, the fraction of background events 
in the normalized-yield difference N„.+ — N„.- estimated 
with Pythia amounts to 20%. The number of ir + events 
after background subtraction is 4510. 

As the recoiling neutron was not detected, t' was de- 
termined from the measurement of the four-momenta of 
the scattered lepton and the produced pion. In order to 
improve the resolution for exclusive events, t' was calcu- 
lated by setting Mj = M n . This offered the possibility 
to discard the lepton energy, the quantity subject to the 
largest uncertainty for the present data sample. This 
procedure results in a factor of two improvement in the 
t' resolution. The same calculation was applied to the 
Pythia Monte Carlo events that were used to subtract 
the background. 

The differential cross section for exclusive ir + produc- 
tion by virtual photons can be written as 



dt' 



1 da ep ^ enn " (x B ,Q 2 ,t' 
^v{xb,Q 2 ) dxBdQ 2 dt't 



(1) 



where Ty is the virtual-photon flux factor. Within the 
Hand convention |30j , this flux factor is 



T v (xb,Q 2 ) 



1 Q 2 l-x B 



Sn M 2 E 2 x 



1 - e 



(2) 



where E is the beam energy and e is the virtual- 
photon polarization parameter. The t' dependence of the 
photon-nucleon differential cross section integrated over 
<j> were extracted from the data as follows: 



dt< ~ Ty ' C Ax b AQ 2 At' k r\ ' 



AjexcL 

JV 7T+ 



(3) 



where N^% cl is the number of n + events after background 
subtraction, k is the probability to detect the scattered 
lepton and the produced ir + within the Hermes spec- 
trometer acceptance, and rj is the radiative correction 
factor. These quantities were determined for each kine- 
matic bin. The symbols Axb, AQ 2 , and At' denote the 
bin size, which was chosen according to the statistical 
precision, instrumental resolution, and kinematic smear- 
ing affecting each individual bin. The t' dependence of 
the differential cross section do ( XB jQ ) was determined 
for four Q 2 bins and the Q 2 dependence of the cross sec- 
tion integrated over t', a(xB,Q 2 ), for three Xb bins. 



The detection probability k was determined using an 
exclusive 7r + Monte Carlo simulation based on a GPD 
model [TBI and a Geant simulation 31 1 of the detector. 
The variable k is the ratio between the number of simu- 
lated events reconstructed in the Hermes spectrometer 
and the number of generated events. These numbers were 
evaluated at the reconstructed and generated kinematics 
respectively and were integrated over the bin size and 
over (j) (and t' for the determination of cf(xb,Q 2 ))- The 
variable k represents the combined effect of the spectrom- 
eter acceptance, the analysis constraints (such as the M\ 
constraint) and the detector efficiencies. For the GPD 
model [l5j| . the power corrections were included and the 
factorized ansatz was used for the t' dependence. How- 
ever, the t' and <f> dependences were modified in order 
to describe the data. Radiative events were included in 
the Monte Carlo simulation according to the code Rad- 
GEN [13] adapted to exclusive meson production using 
the GPD model [HI]. Fig. [3] shows good agreement be- 
tween the kinematic distributions of the data and of the 
exclusive Monte Carlo sample (which is normalized to the 
data by a constant factor) with the exception of yields at 
high values of Q 2 , x and especially— t'. 

The calculated spectrometer acceptance ranges from 
0.1 to 0.7 depending on the kinematic bin, while the effi- 
ciency of the detectors and analysis constraints together 
amounts to about 0.4-0.5. The combination of both leads 
to k values ranging from 0.04 to 0.28. The model depen- 
dence of the determination of k was estimated by using 
the GPD parametrization [HI with different t' and <f> de- 
pendences and by using a different GPD parametriza- 
tion 33j for the Q 2 and xb dependences. It was further 
studied to what degree the relevant kinematic variables 
Q 2 , t' and <f> are correlated by the detector acceptance. 
These model dependences amount to a relative system- 
atic uncertainty in k of less than 15%. 

Events can be smeared from one bin to another by 
multiple scattering and bremsstrahlung. This effect was 
included in the determination of the detection probability 
k. With the selected bin sizes, the fraction of events that 
migrate out of (into) a certain kinematic bin is on average 
12% (15%) and is always below 25% (35%) according to 
the exclusive Monte Carlo simulation. 

The Born cross section was extracted using the radia- 
tive correction factor r], which was determined from the 
ratio of Monte Carlo samples with and without radiative 
effects. The value of r] was found to be 0.77 with little 
variation (less than 3%) as a function of Q 2 , xb or t' with 
the constraints applied to the data to select exclusive 7r + 
production. 

The mean values for xb, Q 2 , and t' for each kinematic 
bin were estimated from distributions generated at Born 
level by the exclusive Monte Carlo simulation. The flux 
factor in Eq. [3] was determined for these mean values of 
Q 2 and Xb- The resulting cross sections were corrected 
for bin-averaging effects to take into account the nonlin- 
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FIG. 3: Distributions of exclusive n + events within the Hermes acceptance as a function of Q 2 , xb, —if and cj> for data (filled 
points) compared with an exclusive Monte Carlo simulation (histogram) based on a GPD parametrization [15I ] which was tuned 
to the data (see text). The data sample is arbitrarily normalized; the Monte Carlo sample is normalized to the data using a 
constant factor. The error bars represent the statistical and systematic uncertainties of the background subtraction added in 
quadrature. 



ear dependence of the cross section within each bin (in 
Eq. [3l a linear dependence inside a bin is implicitly as- 
sumed). These corrections were obtained from the exclu- 
sive Monte Carlo simulation by taking the ratio between 
the cross section evaluated at fixed kinematic values and 
the cross section integrated over the kinematic bins fol- 
lowing Eq. [3l The error arising from the evaluation of 
the flux factor at the mean values is then also corrected. 
The bin-averaging correction factor amounts on average 
to 1.08 and does not exceed 1.2 except for the highest-Q 2 
bin, where it reaches 2.2 at some t' values. 

The integrated luminosity C was determined by com- 
paring the number of inclusive deep-inelastic scattering 
events in the data sample to the yield generated by a 
Monte Carlo sample for inclusive scattering based on 
world data [34|, [35| . The Hermes luminosity detector 
provided another measurement of C, which was used to 
estimate the systematic uncertainty. The integrated lu- 
minosity C amounts to 0.4 fb _1 with a 5% systematic 
uncertainty. 

The total systematic uncertainty of the cross section 
is dominated by the uncertainty of the background sub- 
traction and of the detection probability. The latter takes 
into account the uncertainties due to the model depen- 
dence of its determination and the different detector res- 
olutions in the different data taking periods. 

The measured differential cross section integrated over 
the angle cj) can be written as ^ = ^g- + e^- where 
ot and <tl are, respectively, the contributions of trans- 
versely and longitudinally polarized virtual photons. At 
Hermes the separation of the transverse and longitudi- 
nal components of the cross section is not feasible. How- 
ever, as the transverse contribution is predicted to be 
suppressed by 1/Q 2 with respect to the longitudinal con- 
tribution [l|, the data at larger Q 2 arc expected to be 
dominated by the longitudinal part. Moreover, with the 



presence of the pion pole at low —t', the longitudinal 
part of the cross section is expected to dominate in this 
region. In Fig. 4 and 5, described below, the data are 
compared to calculations for the longitudinal par t of the 
cross section computed using a GPD model hj and to 
calculations for the total cross section computed using a 
Regge model (36|. The cross sections are calculated for 
the mean values of xb, Q 2 -, and t' of the experimental 
data in each bin. 

Fig. [4] shows the t' dependence of the differential cross 
section for four Q 2 bins. As Q 2 is closely related to xb 
(due to the Hermes acceptance as well as the upper 
limit on the pion momentum) low values of Q 2 corre- 
spond to low values of xb- The dashed-dotted lines in 
Fig. [4] show the leading-order calculations of the longi- 
tudinal part computed using the GPD model [lij . The 
GPD E is considered to be dominated by the t— channel 
pion-pole and the pscudoscalar contribution to the cross 
section is parametrized in terms of the pion electromag- 
netic form factor F„. A Regge- inspired t' dependence 
is used for E. The GPD H is neglected here as E is ex- 
pected to dominate the cross section at low —if. The solid 
lines include the power corrections due to the intrinsic 
transverse momentum of the partons and due to the soft- 
overlap contribution, the latter being dominant. While 
the leading-order calculation strongly underestimates the 
data, the calculations including power corrections agree 
with the data for —t' < 0.3 GcV 2 for the four Q 2 bins. As 
the GPD model requires —t' to be much smaller than Q 2 , 
the calculations are not expected to describe the full t' 
range. Furthermore, at larger —if , the data may receive 
a significant contribution from the transverse part of the 
cross section, which is not described by the GPD model. 
Fig. [5] shows the Q 2 dependence of the cross section inte- 
grated over t' for three Xb bins. The Q 2 dependence of 
the data is in general well described by the calculations 
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FIG. 4: Differential cross section for exclusive 7r + production by virtual photons as a function of —t' for four Q 2 bins. The 
inner error bars represent the statistical uncertainties and the outer error bars represent the quadratic sum of statistical and 
systematic uncertainties. The curves represent calculations based on a GPD model [15j for ^f- using a Regge-type ansatz for 
the t' dependence (dashed-dotted lines: leading-order calculations, solid lines: with power corrections) and a Regge model [36l ] 
for 4f (dashed lines) and (dotted lines). 



A 0.02 <x B < 0.15 
• 0.15 <x B < 0.26 
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FIG. 5: Cross section for exclusive % + production by vir- 
tual photons as a function of Q 2 for three xb bins and inte- 
grated over t' (colour online). The inner error bars represent 
the statistical uncertainties and the outer error bars repre- 
sent the quadratic sum of statistical and systematic uncer- 
tainties. The curves represent calculations based on a GPD 
model [l5j for ol using a Regge-type ansatz for the t' depen- 
dence (dashed-dotted lines: leading-order calculations, solid 
lines: with power corrections) and a Regge model [3(| for a 
(dashed lines). The thin, medium and thick lines correspond 
to the low, medium and high xb values, respectively. 



from the GPD model [15[ with the inclusion of the power 
corrections, although the magnitude of the cross section 
is underestimated. The data support the order of magni- 
tude of the power corrections for the calculations by the 
GPD model [l5| at low — t', a region where the longitu- 
dinal part of the cross section is expected to dominate, 
and for the available Q 2 range. 

Both the transverse and longitudinal parts of the cross 
section were computed using a Regge model (36j . where 



pion production is described by the exchange of tt and p 
Regge trajectories. In this formalism, the meson-nucleon 
coupling constants are fixed by pion photoproduction 
data. In the original version [20|], the irir~f form factor 
is fixed by pion form factor measurements, while the 71737 
transition form factor is unconstrained. In the version 
used here (3(|, both form factors are taken to be both 
Q 2 - and ^'-dependent. The dashed (dotted) lines in Fig. [4] 
show the total (longitudinal) cross section computed us- 
ing the Regge model. In this model the transverse part of 
the cross section is estimated to represent from 6% to 8% 
of the total cross section at — t' = 0.07 GeV 2 and from 
15% to 25% of the total cross section integrated over t', 
confirming the expected suppression of the transverse to 
the longitudinal part of the cross section. However data 
from JLab 25j at lower center of mass energy (W 2 = 4.9 
GeV 2 ) show that the transverse part of the cross section 
is underestimated by the Regge model by a factor of three 
to four. It is not clear if this also holds at the higher 
W 2 and —t' values of the Hermes data. Compared to 
the calculations for the longitudinal cross section from 
the Regge model (dotted lines) , the t' dependence of the 
GPD model [IH (solid lines) in Fig. |4] appears too steep. 
The total cross section computed by the Regge model 
describes well the t' dependence of the differential cross 
section (dashed line on Fig. [3]) and the Q 2 dependence 
of the cross section integrated over t' (dashed lines on 
Fig. El). The model calculations give also a good descrip- 
tion of the magnitude of the data except at low —t' for 
Q 2 < 3 GeV 2 , where the calculations overestimate the 
data up to 70%. 

In conclusion, the cross section was measured for ex- 
clusive electroproduction of n + mesons from a hydrogen 
target as a function of —t' for four Q 2 bins and as a func- 
tion of Q 2 for three xb bins. A model calculation for 
the longitudinal part of the virtual-photon cross section 
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based on the Generalized Parton Distributions [15[ does 
not describe the data at the leading order in the present 
Q 2 range. However, if power corrections are included, 
the model is in fair agreement with the magnitude of 
the data at low values of — t'. In this kinematic region, 
where ir + production by longitudinal photons is expected 
to dominate, the data support the order of magnitude 
of the power corrections. A model calculation based on 
the Regge formalism for both the l ong itudinal and the 
transverse part of the cross section [361 ] provides a good 
description of the magnitude and the t' and the Q 2 de- 
pendences of the data. 
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gratefully acknowledge the Desy management for its 
support, the staff at Desy and the collaborating institu- 
tions for their significant effort, and our national funding 
agencies and the EU RII3-CT-2004-506078 program for 
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